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Kelvin waves (kelvons) — helical waves on quantized vortex lines — are the normal modes of vortices 
in a superfluid. At zero temperature, the only dissipative channel of vortex dynamics is phonon 
emission. Starting with the hydrodynamic action, we derive the Hamiltonian of vortex-phonon 
interaction, thereby reducing the problem of the interaction of Kelvin waves with sound to inelastic 
elementary excitation scattering. On the basis of this formalism, we calculate the rate of sound 
radiation by superfluid turbulence at zero temperature and estimate the value of short-wavelength 
cutoff of the turbulence spectrum. 
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Since the early study of phonon scattering by vortex 
lines much interest has been attracted by the prob- 
lems of the interaction of Kelvin waves with density- 
distortion modes @,0,0,IHE10]- In helium, Kelvin waves 
are a fundamental ingredient of the evolution of super- 
fluid turbulence H \£M H G3 El A Kolmogorov- 
type cascade of Kelvin waves has been argued 0, fiol 
to be responsible for the decay of superfluid turbulence 
at zero temperature ^3 • As proposed by Vinen 0, , 
the sound radiation by short-wavelength kelvons leads to 
a dissipative cutoff of the Kolmogorov energy flux. In 
neutron stars, the excitation of Kelvin waves due to the 
interaction with the nuclei in the solid crust is suggested 
to be the main mechanism of pulsar glitches 2\. Not 
long ago, quantized vortices were created in atomic Bosc- 
Einstein condensate (BEC) 13j. Since then, nonlinear 
kelvon dynamics have become an attractive topic in the 
field of ultra-cold gases 0, IE lljl - Recent experiment 
in gaseous BEC by Bretin et. al. 5] and the subsequent 
numerical simulation by Mizushima et. al. describe 
the excitation of Kelvin waves by the quadrupole modes. 
In their theoretical study, Martikainen and Stoof Q ob- 
tained the interaction Hamiltonian of kelvons with the 
quadrupole modes in a model of a vortex line in a stack 
of two-dimensional BEC. 

We develop a systematic approach to the problem of 
interaction of phonons with vortices in the hydrodynamic 
regime, i.e. when any physical length scale is much larger 
than the vortex core size clq, which allows one to describe 
vortices as geometrical lines $] . We derive the interaction 
Hamiltonian basing the analysis on the small parameter 
/3 = aak -C 1, where k is the largest wave number among 
kelvons and phonons. To employ the transparent descrip- 
tion in terms of the normal modes, we confine ourselves 
to the case of weak nonlinearity. For kelvons this implies 
that the amplitudes &/. of the Kelvin waves of the typical 
wavelength A ~ k^ 1 are much smaller than A, which is 
expressed by the small parameter — b^k <C 1. For 
phonons this requires that r\ -C n, where r\ is the number 
density fluctuation in a sound wave and n is the average 
number density. The obtained result allows us to rigor- 
ously describe the radiation of sound by kelvons, which 
we apply to the problem of superfluid turbulence decay at 



zero temperature. We should mention that, in superfluid 
turbulence, kelvons are actually superimposed on vortex 
kinks with curvature radii much larger than kelvon wave- 
lengths We shall neglect such large-scale curvature 
and address this issue in the end of the paper. 

The condition a$k ~ (3 <C 1 implies that the typical 
vortex line velocities are much smaller than the speed 
of sound. Along with rj <C n it leads to the fact that 
the vortex-phonon coupling contributes only small cor- 
rections to the dynamics of the non-interacting vortex 
and phonon subsystems. Therefore, a perturbative ap- 
proach is applicable, provided the interaction energy is 
written in terms of the canonical variables. Normally, 
the form of the canonical variables comes from the so- 
lution of the interaction-free dynamics. However, when 
studying vortices separate from phonons,one naturally 
neglects the compressibility of the fluid [g, since finite 
compressibility leads only to higher-order "relativistic" 
corrections. As a result, the dynamics of vortices are de- 
scribed by the Hamiltonian, written in terms of the geo- 
metrical configuration of the vortex lines ^3 • When the 
vortices are absent, the phonon modes come from the bi- 
linear Hamiltonian for the density fluctuation 77(1*, t) and 
the phase field (p(r,t), which determines the velocity in 
the density wave (see, e.g., HH). If finite compressibil- 
ity of a superfluid is taken into account in order to join 
the subsystems, the positions of vortices and the fields 
?7(r, t), <p(r, t) are no longer the sets of canonical variables 
because of the variable-mixing term in the Lagrangian. 
Introducing the interaction, one has to simultaneously 
reconsider the canonical variables. 

The small parameters allow us to obtain an asymp- 
totic expansion of the canonical variables by means of a 
systematic iterative procedure. Physically, the procedure 
restores the retardation in the adjustment of the super- 
fluid velocity field to the evolving vortex configuration. 
It qualitatively changes the structure of the Hamiltonian 
with respect to the terms responsible for the radiation of 
sound and the relativistic corrections to the vortex dy- 



Long-wave superfluid dynamics at zero temperature 
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are described by the Popov's hydrodynamic action [l6| 
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dtd 3 r 
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(1) 



Here the spatial integral is taken over the macroscopic 
fluid volume, $(r, t) is the phase field, which determines 
the velocity according to v = (1/too)V$ (ft = 1), mo 
is the mass of a particle, >c is the compressibility, the 
dot denotes the derivative with respect to time, and the 
vortex core radius is given by ao ~ yjxjnm^. 

The phase $ is non-single-valued and contains topo- 
logical defects, the vortex lines. The velocity circulation 
around each vortex line is quantized: 



V$(r) • dv = 2tt. 



(2) 



The defects can be separated from the regular contribu- 
tion: 



$ = $ +ip, 



(3) 



where $o is non-single- valued and satisfies ? while ip is 
regular and \/ip is circulation-free. The standard decom- 
position into vortices and phonons is done by introducing 
an additional constraint that to the zeroth approximation 
eliminates the coupling between $o and ip in the Hamil- 
tonian, namely 



A$o(r) =0. 



(4) 



(Physically, this parametrization is suggested by the 
velocity potential of an incompressible fluid.) With 
Eqs. ©,0 the Lagrangian becomes 



d 3 i 



-n4> — rjcp — ry$ 



H, 



where H = H v . 
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d 3 r V$ 



H ph = d 3 r 



H ^ ~ 2m 



d A r 
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j) V$ + 2r?Vv? • V$ 



(5) 

(6) 
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The coupling between the vortex variable, <E»o : and the 
density waves, {i],ip}, is determined by H[ ai and the 
time derivative term Jd 3 r f^Oj both being first-order cor- 
rections to the non-interacting parts. Following stan- 
dard perturbative procedure, we first neglect this cou- 
pling to find the non-interacting normal modes. The 
vortex part of the Lagrangian is then given by L vor = 



;jd 3 r $ 



For the sake of simplicity, from 



now on we consider a solitary vortex line; the generaliza- 
tion is straightforward. Let the two-dimensional vector 
Po( z o) = { x o( z o), Vo{ z o)i 0) describe the position of the 
vortex line in the plane z = zq of a Cartesian coordinate 



system, where the ^-direction is chosen along the vortex 
line. The field $o is a functional of Po(z), hence 



d 3 r <F 



o 



dz po ■ — jd 3 r $ . 



(9) 



To obtain J d 3 r $ , it is sufficient to calculate the integral 
/ d 3 r 5&o, where <5<I>o is the variation of the phase field 
due to the distortion of the vortex line by Spo(z). Using 
the identity A(p 2 ) = 4, with p = (x,y,0), and Eq. |@J, 
obtain 



Jd 3 r 6% = i Jd 3 r V • [ <5<f V(p 2 ) ] 



(10) 



The variation 5<f>o(r) can be viewed as being produced by 
two vortex lines with opposite circulation quanta and sep- 
arated by Spo(z). In view of Eq. (O the field <5$o( r ) ex- 
periences a jump of 2ir across the surface S that extends 
between these vortex lines along the vector field 5p a (z), 
therefore the integration volume must have a cut along 
S. Applying the Gauss theorem, we rewrite Eq. (|lTj|) as 
the surface integral (1/2) <f s <5$ (p • dS) yielding 



d 3 r $ 



dz [z x po(z)] ■ po(z) 



(11) 



Introducing the complex variable w(z) = 
\J nmon/2 \x{z) + iy(z)], where k = 27t/too is the 
velocity circulation quantum, obtain 



\ j dz [' 



iw w — IW w 



\w, w 



(12) 



which implies that w(z) and w*(z) are the canonical vari- 
ables with respect to L vor . The energy rewritten in 
terms of w{z) and w*(z) gives the vortex Hamiltonian 
|10| . In this paper we keep only the bilinear term of the 
expanded with respect to < 1 Hamiltonian, the re- 
maining terms determine the interaction of kelvons |l2j | . 
The result is £3 
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e k = ( K /4^)ln(l/fca )fc 2 , (13) 



where and a\. are the kelvon creation and annihilation 
operators, obtained, assuming the periodicity along z, by 
the Fourier transforms = L^ 1 / 2 f w(z) exp(-ikz) dz, 
where L is the system size in the ^-direction, and w(z) is 
understood as a quantum field. 

The sound waves are described by the Lagrangian 
£ P h = Jd 3 r(-r)tp) - Hphiv,^], with H ph given by Q. 
We assume that the system is contained in a cylinder 
of radius R with the symmetry axis along the z-direction 
and that the system is periodic along z with the period L. 
In the cylindrical geometry, the phonon fields r](r,9,z), 
ip(r, 9, z) are parametrized by phonon creation and anni- 
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hilation operators c s , c\ as 

V = Y^ [ Xs c s + x* s 4 ] ) 

p = -i^ *JTj2u^€ [ Xs c s - Xt 4 ] , (14) 

S 

Xs = Xs(r, 0, z) = K mqr {r) Y m (6) Z q ,(z), 

where TZ mqr (r) = (nq r /R) 1 / 2 J m (q r r), Y m (6) = 
(27r) -1 / 2 exp(im^), 2 qz (z) = L~ x l 2 exp(iq z z), s stands 
for {q ri m, q z }, and J m (x) are the Bessel functions of the 
first kind. The phonon Hamiltonian then reads 



(15) 



with q = y/ q 2 + g 2 anc [ c — ^/ n/ xttlq. 

Now we address the coupling between phonons and 
vortices. In terms of the obtained variables, the La- 
grangian ||SJ takes on the form 



L 



E 



i d k a\ 



i c s c\ — T — H, 



(16) 



where T = Jd 3 r r/$ = T{a k ,a kl a\,a\, c s ,c s ,c\,c\}, 
and H — iJ vor + H ph + H[ nt . The coupling term T plays 
a special role in the Lagrangian (|16[) . This term is linear 
in time derivatives of the variables and thus can not con- 
tribute to the energy in accordance with the Lagrangian 
formalism. Moreover, because of the time derivatives in 
T, the equations of motion in terms of {a k ,a k },{c s ,cl} 
take on a non-Hamiltonian form. This implies that the 
chosen variables become non-canonical in the presence 
of the interaction, and therefore the total energy, in 
these variables can not be identified with the Hamilto- 
nian. 

There must exist such a variable transformation 
{a k , <4}, {c s , c|} — > {d k , a\}, {c s , c|} that restores the 
canonical form of the Lagrangian, L 



ic s c\ 



H{d k , flfc, Cs, c\}. The canonical variables are 



obtained by the following iterative procedure. The term 
T is expanded with respect to at <C 1, /3 <C 1 and 
7] < n yielding T = TW + T< 2 ) + ■ • • . Then the vari- 



aji , c., 



ables are adjusted by a k — * a k 

where a^\{a k ,a\,c s ,c\}) and d 1 ' ({a k , a\, c s , c\}) are 
chosen to eliminate the term T^ 1 ' in j!6|>. As a result, 
T — > 0+T'( 2 ) + - • • , where the prime means that the struc- 
ture of the remaining terms has changed. At the next 
step, T ,{ - 2 > is eliminated by a k — > a k + a k , c s — > c s + c s 
and so on. By construction, the canonical variables are 

u ~ , (!) , (2) , - i (i) , (2) , 
given by a k = a k +a k +a k ' -\ , c s = c s +c s '+c s M , 

and likewise for the conjugates. In practice, only the first 
few terms are enough, as the rest ones give just higher- 
order corrections. 

The explicit expression for T is obtained following the 
steps of the derivation of Eq. (|11|) . The only difference 



here is that the role of the auxiliary function p 2 is played 
by Q, defined by AQ(r) = -q(r): 



T = 2ir dz [z x VQ(p (z), z)] ■ p Q (z). 



(17) 



In accordance with q r po ~ ftkpo <C 1, we expand 
the radial functions, retaining only the greatest term for 
each particular angular momentum to. Switching to the 
cylindrical coordinates, po = (po cos 7, po sin 7, 0), and 
noticing that 



7Po 



dr p d9 



H mqr (r)Y m (9) 



r=p ,tf=7 

d{w m )/dt, m>0 
d(w*^)/dt, m < 



(18) 



obtain 



...k m C s -jMk^ 



• a k m ) 



s,k±...k 

a . + 

- iB s ^ kl ... km c s -^i a ki •••"■k 
where the sum is over all s with to ^ and 
.4 



H.c. 



(19) 



As = 



s,kl...k m — — ©(Wl) A s 5 kl + ...+ km .-qz-, 

B sM ... km = (-l)HG(-m- 1) A s S kl+ ... +km ^, 

, i-M , , 2-|m| \ n 

a q 



2-Jm] | m |+i _ 2 

n 2 (wo K J 2 9r 9 



where 0(m) 



2|m|/2+l | m |J ^(|m|-l)/2 JJ1/2 

1, m > 



, (20) 



0, TO < 



Thus, 



a k — a k 



(21) 



(we omitted the terms that do not contain phonon op- 
erators and thus result only in relativistic corrections to 
the kelvon spectrum and kclvon-kelvon interactions), 

c s = Cs + (1 - #ro,o) E [^s.fei-fem^li ■ ■ ■ "fe m 



+ B s , kl ... km d kl . . . d kn 



s = {q r ,m,q z }. (22) 



The Hamiltonian H is then given by the energy (JBJ- 
(JSJ) in terms of the variables {a k ,a k },{c s ,cl}. Up to 
neglected relativistic corrections, the variable transfor- 
mation does not change the spectrum of the elementary 
modes: the zero-order Hamiltonians are given by l|13f) . 
(TT5)) in terms of {d k ,d\}, {c s ,c|}. The transform 
applied to (|15f) generates the interaction term 



H 



(rad) 



(1 — ^ ro ,o) cu s A St h 1 ...k m o/ ki 

s,{ki} 



~t ~t 



+ u s Bs :kl ... km d kl 



■ak m c\ 



+ H.c. (23) 
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Remarkably, the energy term oc Jd 3 r rj\ V$o | in JSJl, 
which results in the same operator structure as Eq. (|23[1 . 
is irrelevant, being smaller in (3 <C 1. It can be checked 
straightforwardly, that the term oc Jd 3 r nVip ■ V$o m 
Eq. JSJl gives Fetter's amplitudes of the elastic and inelas- 
tic scattering of phonons (see, however, the remark at 
the end of the paper). In addition, this term leads to a 
macroscopically small splitting of the phonon spectrum 
due to the superimposed fluid circulation. 

A kelvon is known to carry a quantum of (negative) 
angular momentum projection 0]. The interaction H23|) 
explicitly conserves the angular momentum: a real pro- 
cess of the emission of a phonon with the angular mo- 
mentum (— m) requires an annihilation of m kelvons. 
Since e k ~ (aok)uj k , the total momentum transferred 
to phonons in a radiation event should be small in or- 
der to satisfy the energy conservation. Thus, the radi- 
ation by one kelvon is kinematically suppressed. The 
leading radiation process is the emission of the m = — 2 
(quadrupole) phonon mode, the events involving more 
than two kelvons being suppressed by a k <C 1. First- 
order processes of the two-phonon emission come from 
the term oc Jd 3 r nVip ■ V^o of ijSJl. The amplitude of 
these processes is suppressed by the relativistic parame- 
ter (3 < 1. 

The Hamiltonian l)23|) can be used to obtain the rate 
of sound radiation by superfluid turbulence at zero tem- 
perature. At large wave numbers, where the radia- 
tion is appreciable, Kelvin- wave turbulence is charac- 
terized by a k <C 1 and features the spectrum n k = 
{a\,a k ) oc fc~ 17 / 5 [l^. The occupation number decay rate 
is h k = - J2sM Ws.fe.fei, where W s ,u,ki is the probability 
of the event |0 S , n k , n kl ) — > |l s , n k -l, n kl -l) per unit 
time. Applying the Fermi Golden Rule to W Sykykl with 
the interaction l)23|) and replacing the sums by integrals, 
obtain 



(moK/27r) E 
157rnmo 



\n 5 (l/a k)(a k) 5 k 5 n 2 k . (24) 



Following Ref. |TJ|, this formula allows us to obtain the 
momentum scale £; p h, at which the kelvon cascade is cut 
off by the radiation of sound: 



« P h 



[a /R 



6/31 



[\n(R /a )] 



24/31 



"•■0 ' 



(25) 



where Rq is the typical distance between the vortex lines 
in the tangle. 



Finally, we comment on Vinen's estimate of the power 
radiated per unit length of the vortex line, Eq. (2.24) 
of Ref. 4]. The power at the momentum k can be de- 
termined by IK = -J2k>~ k £ k>h k >/L, with n k > given by 
l)24p. In Ref. [4(, the retarded potential method is em- 
ployed giving the estimate Yi' k oc b\ oc n k , while Ilfc oc n\. 
In the quasi-particle language, 11^ implies that the radi- 
ation is governed by the conversion of one kelvon into a 
phonon. Vinen argues that this process becomes allowed 
in superfluid turbulence, where kelvons are actually su- 
perimposed on vortex kinks of typical size ~ Rq 3> k , 
or, equivalently, the kelvon coupling to a kink lifts the 
ban on single-kelvon radiative processes by effectively re- 
moving the momentum conservation constraint. We note 
that the probabilities of such elementary events are likely 
to be suppressed exponentially, as it is gencrically the 
case, say, for soliton-phonon interactions (see, e.g., 01 
and references therein.) The processes involving kelvon- 
kink coupling should contain an exponentially small fac- 
tor ~ exp(— i?ofc), which arises from the convolution of 
the smooth kink profile with the oscillating kelvon mode. 



In the problem of phonon scattering from a vortex, 
the answer depends on whether the vortex is pinned or 
free [l^. In fact, the elastic scattering matrix element 
of Ref. corresponds to the pinned case only. In the 
Hamiltonian formalism developed here, the difference is 
due to an additional matrix element generated by the 
|m| = 1 term of the interaction Hamiltonian (|23[) in the 
second order of perturbation theory. We are indebted to 
Edouard Sonin for raising this question. Furthermore, 
a second-order amplitude generated by a combination of 
the \m\ = 2 term with the \m\ = 1 one is of the same or- 
der as the direct first-order amplitude of inelastic phonon 
scattering. Thus, the result of Ref. for inelastic scat- 
tering needs to be corrected accordingly. 
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